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An acoustic transfer function relating combustion noise and turbine exit noise in the presence of 
enclosed ambient core noise is investigated using a dynamic system model and an acoustic 
system model for the particular turbofan engine studied and for a range of operating conditions. 
Measurements of cross-spectra magnitude and phase between the combustor and turbine exit and 
auto-spectra at the turbine exit and combustor are used to show the presence of indirect and direct 
combustion noise over the frequency range of CM-00 Hz. The procedure used evaluates the ratio of 
direct to indirect combustion noise. The procedure used also evaluates the post-combustion 
residence time in the combustor which is a factor in the formation of thermal NO v and soot in this 
region. These measurements are masked by the ambient core noise sound field in this frequency 
range which is observable since the transducers are situated within an acoustic wavelength of one 
another. An ambient core noise field model based on one and two dimensional spatial correlation 
functions is used to replicate the spatially correlated response of the pair of transducers. The spatial 
correlation function increases measured attenuation due to destructive interference and masks the 
true attenuation of the turbine. 
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NOMENCLATURE 

B e Effective resolution bandwidth, B e = l/T d = rJN, 
2 Hz 

C Cost function 

f u Upper frequency limit, /„= \/2At = rJ2, Hz (32 
768 Hz) 

G x Jf) Estimated input signal one-sided auto spectral 
density function 

G xy (f) Estimated one-sided cross-spectral density 
function 

Gy t y(f) Estimated output signal one-sided autospectral 
density function 

H(f) Frequency response (transfer) function 

K Turbine attenuation or producton/creation factor, 

£2 

N Segment length, number of samples per segment, 
N=r s T d , (32 768) 

N d Direct combustion noise signal system noise 
Nj Indirect combustion noise signal system noise 
N t Tailpipe signal system noise 
n a Number of segments/blocks with 50 percent over- 
lap, n 0 = 2 B e T total « 273 

n s Number of disjoint (independent) segments used 
in spectra estimates, n s = B e T t ota i ~ 128 
r m Microphone separation 
r s Sample rate, samples/s (65 536) 

S Cost function sums 

T total Total record length, T tota ] = n s T d , s («70 s) 

T d Record length of segment, T d = N/r s , 0.5 s 
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U x, y (f) Aligned cross-spectral density function 

x{t) Input signal from combustor pressure sensor CIP1 
y(t) Output signal from turbine exit pressure sensor 

10 Turbine exit signal, T551 

1 1 Turbine exit signal, T552 

9 Combustor probe signal, CIP1 
c Core 

d Direct combustion noise 
e Effective 

i Indirect combustion noise 
o Zero frequency 
pc Post-combustion 
TE Turbine exit 

x Input signal 

y Output signal 

oc c i,a c2 One dimensional core noise spatial correlation 
coefficient level constant, eg,en 
o. TE Two dimensional turbine exit noise spatial corre- 
lation level constant, e 6 

1 1 Angle of incident plane wave to a pair of sensors 

At Sampling interval, At = 1 /r s s 
e Model parameter 

)/ Two dimensional turbine exit noise spatial corre- 
lation coefficient function 

v Ratio of contribution of direct combustion noise 
to indirect combustion noise, e\ 
p(f) Spatial correlation for two separated turbofan 
engine core pressure sensors 
t 0 Indirect combustion noise time delay alignment, 
3.05 ms 

t a Acoustic wave travel time between sensors, ms 
t c i,t C 2 One dimensional core noise spatial correlation 
coefficient time constant, € 9 , 612 , ms 
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t pc Post-combustion residence time, 64, ms 

x TE Two dimensional turbine exit noise spatial corre- 
lation time constant, 67, ms 
x T Turbine frequency response time constant, e 3, s 
( p Post-combustion residence time phase shift, deg 

cp 0 Zero frequency post-combustion residence time 
phase shift, e 5 , deg 

i 9 Cl , i) Cl Zero frequency one dimensional core noise spatial 
correlation coefficient phase shift, eio,ei3, deg 
£1,1^2 One dimensional core noise spatial correlation 
coefficient function 
A Acoustic 
E Estimated 
M measured 

I. INTRODUCTION 

Low frequency noise generated in the turbofan engine 
core by direct and indirect combustion noise may make a 
significant contribution to the overall noise signature in the 
aft direction at the low power settings used on an airport 
flight trajectory. In addition, low frequency core noise may 
also be a problem for future aircraft using lean-burning 
ultra-low-emissions combustors. The results of using pres- 
sure sensor measurements in a turbofan engine combustor 
and at the turbine exit to determine a low frequency acoustic 
linear system transfer function relating the turbine exit noise 
to the combustion noise are reported here. The linear system 
transfer function model is enhanced by a spatial acoustic cor- 
relation function between the sensors in order to account for 
the enclosure sound field. Because proprietary geometry and 
flow information is not available, the method used avoids the 
usual limitation imposed by the need to know the exact ge- 
ometry and flow conditions in the engine core to account for 
the acoustics of the enclosure. The approach used is based 
on an approximate analytical formulation. In the frequency 
range studied, indirect combustion noise is separated from 
direct combustion noise since it has a phase relationship to 
direct combustion noise characterized by a post-combustion 
(post-flame) residence time delay. 

Two possible low frequency turbofan engine noise sour- 
ces are “direct” and “indirect” combustion noise. The source 
of combustion noise attributed to the unsteady pressures pro- 
duced by the unsteady combustion process that propagate at 
the speed of sound through the turbine to the far field is called 
the “direct” combustion noise source. The generation of 
direct combustion noise from a open non-premixed turbulent 
flame is discussed by Ihme and Pitsch. 1 The generation of 
direct combustion noise in a combustor is discussed by Liu 
et al? The other source of turbofan engine combustion noise 
is known as the “indirect” mechanism in which the noise is 
generated in the turbine by the interaction of entropy fluctua- 
tions, which also originate from the unsteady combustion 
process but propagate at the flow velocity. Noise is generated 
as entropy fluctuations propagate through regions character- 
ized by mean flow velocity gradients or pressure gradients in 
the turbine stages. A summary of early references using 
“actuator disc” theory to study turbine generated indirect 
combustion noise is given by Cumpsty. Bodony presents 
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computational and analytical results on the interaction of a 
convected entropy disturbance with an air-foil shaped body. 4 
A related topic is the noise generated by propagation of an 
entropy disturbance through a nozzle considered theoretically 
by Marble and Candel. Leyko et al ., 6 and Howe. 7 Experi- 
mental studies of indirect combustion noise generated by 
flow of an entropy disturbance through a nozzle have been 
conducted by Bake et al . 8 and Rausch et al . 9 A theoretical 
study of the generation of indirect combustion noise in an 
auxiliary power unit (APU) has been conducted by Xu, Tam, 
and Parrish. 10 Combustion chamber instability can be caused 
by a feedback mechanism which couples combustion cham- 
ber acoustics with convectively transported fluctuations of 
entropy as noted by Polifke, Paschereit and Dobbeling. 1 1 

Investigations of indirect and direct turbofan engine 
noise in the far field were conducted by turbofan engine 
companies with NASA sponsorship. Experimental studies of 
indirect combustion noise from a Pratt & Whitney PW 4098 
(East Hartford, CT) turbofan engine using an aligned and 
purposely nonaligned coherence procedure on acoustic 
measurements in the engine combustor and in the far field 
are discussed by Miles. 12 ' 11 Experimental studies of indirect 
combustion noise from a Honeywell TECH 977 (Phoenix, 
AZ) turbofan engine using cross-spectra and correlation 
function procedures to estimate the post-combustion (post- 
flame) residence time from acoustic measurements in the 
engine combustor and in the far field are discussed by 
Miles. 14-1 7 The post-combustion (post-flame) residence time 
is the time taken for an entropy disturbance to travel in the 
combustor from the flame to the turbine entrance at the flow 
velocity. These studies showed indirect combustion noise 
was present in these turbofan engines at low frequencies. 

The relative importance of direct and indirect combus- 
tion noise mechanisms in creating core noise in modern 
engines and future engines is unknown. In addition, the fun- 
damental combustion noise prediction capability of the 
NASA Aircraft Noise Prediction Program (ANOPP) is based 
on physical models attributing the core noise to direct com- 
bustion noise. In order to assess and improve the NASA 
combustion noise prediction capability, noise source separa- 
tion studies using pressure sensors in the combustor and at 
the turbine exit have been made by Hultgren and Miles. 18 
Hultgren 19 used noise separation signal processing techni- 
ques with three engine-internal sensors to determine the 
turbine transfer function. Hultgren 19 also suggested an 
improvement for the direct combustion-noise module in 
ANOPP. 

A lot of research on topics related to direct and indirect 
combustion noise has been performed without examination 
of turbofan engine core acoustic measurements to study 
direct and indirect combustion noise as generated in a turbo- 
fan engine. As a step toward improving the understanding of 
the acoustics in the combustor-turbine system, an acoustic 
linear system transfer function model relating direct and 
indirect combustion noise and the turbine exit noise in the 
presence of ambient core noise is reported here for a range 
of operating conditions. Using two engine-internal sensors 
from the same test program studied by Hultgren, 19 the 
combustor and turbine exit auto-spectra and the combustor/ 
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turbine cross-spectra are calculated and are used in develop- 
ing the transfer function model reported here. The original 
model was based on an acoustic linear system transfer func- 
tion between the combustion noise and the turbine exit noise 
and did not take into account the effect of the enclosure. The 
result was a noticeable gap between the measurements and 
the model results. The advance demonstrated here is that by 
providing an enclosure sound field based on spatial correla- 
tion in the engine core using an approximate analytical for- 
mulation, the resulting model is in improved agreement with 
measured acoustic cross-spectra phase. The results presented 
by Hultgren 19 do not account for an ambient enclosure sound 
field based on spatial correlation in the engine core. This pa- 
per shows that by including an ambient enclosure sound field 
the large scale distortion observed in the cross-spectrum 
phase angle can be replicated. In addition, by including an 
ambient enclosure sound field the reduction in the measured 
auto and cross-spectrum amplitude at low frequencies is also 
replicated. These results are of minor consequence to the 
conclusions drawn by Hultgren. 19 However, they radically 
complicate any simple approach one might propose to inves- 
tigate the relation between combustion noise and noise at the 
turbine exit that is based on comparing calculated results 
with measurements. 

In order to gain physical insight into combustion noise, 
auto and cross spectra measurements made using a pressure 
sensor in the combustor and pressure sensors at the turbine 
exit have been studied to determine how much noise at the 
turbine exit is due to indirect combustion noise and how 
much noise is due to direct combustion noise. In addition, 
the cross-spectrum phase angle slope can provide informa- 
tion on the post-combustion residence time in the combustor 
which is a factor in the formation of thermal NO, and soot in 
this region. Unfortunately, the low frequency measurements 
are made in a confined space and the measured auto spectra 
at the turbine exit and the cross-spectra are masked by an in- 
terference pattern due to the ambient sound field which can 
be measured because the pressure sensors are situated within 
an acoustic wavelength of one another. The core noise pair 
of sensors are within an acoustic wavelength of one another 
and behave in a manner similar to a pair of sensors in the 
ocean or in a diffuse sound field in a reverberant room or 
wind tunnel when a pair of sensors in these venues are within 
an acoustic wavelength of one another. 20-26 One dimen- 
sional and two dimensional spatial correlation function mod- 
els are used to study the correlated responses of the 
transducers to the ambient acoustic core noise field. The de- 
velopment of a simple physical model for the transfer func- 
tion lead to the discovery of spatial correlation in the 
enclosed ambient core noise field. 

Turbine noise was studied extensively in the 1970s. 
Studies were conducted on the topic of turbine tone genera- 
tion. Other studies deal with the topic of attenuation of direct 
combustion noise. The topic of indirect combustion noise 
was studied analytically and in model scale tests. Indirect 
combustion noise was not thought to be an important con- 
tributor to turbofan engine core noise. Consequently, experi- 
mental studies of indirect combustion noise in turbofan 
engines were not conducted in this time period. Presently, as 
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a result of previous studies, core noise from turbofan engines 
is treated as being due to such individual sources as the com- 
pressor, combustor, and turbine. The turbine in previous ana- 
lytical and experimental studies has been treated as a device 
that attenuates the direct combustion noise with an attenua- 
tion constant that depends only on the turbofan engine opera- 
tion parameters such as the pressure and temperature at the 
inlet and exit of the turbine and not on frequency. Stone, 
Krejsa, and Clark 27 discuss the current status core noise 
modeling for turbofan engines and provide a historical 
perspective. 

The net travel time of the indirect combustion noise sig- 
nal from the combustor to the far field is increased since the 
travel velocity of the entropy fluctuations to the turbine is 
the flow velocity in the combustor. This flow velocity is a 
small fraction of the speed of sound. Miles el al. ' J ‘ have 
shown that the pressure and entropy should be in phase in 
the combustor. Consequently, one might expect that the 
pressure signal from an indirect combustion noise source 
would be delayed relative to a pressure signal from a direct 
combustion noise source since an indirect combustion noise 
signal does not travel with the speed of an acoustic wave 
until it interacts with the turbine. Miles 14-16 with data from 
the Honeywell TECH 977 engine test program shows that 
the cross-spectra and correlation function between a combus- 
tor sensor and far-field microphones are tools that provide a 
way to separate “direct” and “indirect” coherent combustion 
noise due to this travel delay time. This paper uses measure- 
ments in the combustor and at the turbine exit from the Hon- 
eywell TECH 977 engine test program. This paper uses tools 
that are part of signal-processing theory to study a combustor 
pressure sensor signal and a turbine exit pressure sensor sig- 
nal. The cross-spectral density phase measurement identifies 
a time delay that corresponds to the convective time delay 
which identifies the presence of indirect combustion noise. 
This time delay is the post-combustion residence time in the 
combustor which is a factor in the formation of thermal N Ch- 
an d soot in the post-combustion zone in the combustor. The 
coherence between the two sensors identifies the spectral 
region of importance to indirect combustion noise to being 
in the 0-400 Hz frequency range. This paper discusses the 
cross-spectral density phase angle and the coherence. 

The investigation in this paper covers the frequency 
range below 400 Hz and not the turbine tone source region 
in the range of 3 to 9 kHz which has been studied experi- 
mentally by Enghardt et al. 30 and computationally by Van 
Zante and Envia. 3 In addition, compressor noise is above 
the 400 Hz frequency range 32 and does not contribute in the 
frequency range considered. 

This source identification problem is similar to the pa- 
rameter or system identification problem. Most recent stand- 
ard treatments of system identification problems use the time 
domain and lump together linear and nonlinear system treat- 
ments. These modern methods can account for changes in 
dynamics as a function of excitation amplitude. Unfortu- 
nately, procedures discussed assume one has control of the 
system input signals so one can study system response to dif- 
ferent amplitudes. This ability is not generally available for 
a turbofan engine noise source study because abrupt changes 
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or perturbations introduced might for example modify the 
pattern factor which describes the spatial temperature distri- 
bution across the turbine blades. This can cause enough 
stress on rapidly rotating blades to create turbine blade stress 
fractures and destroy the engine. In addition, these modern 
methods do not cover systems with time delays, such as the 
convective time delay connected with the indirect combus- 
tion noise. As an alternative, one approach is to create a lin- 
ear model approximation that can be considered to be the 
first term in a series of functions that would include nonlin- 
ear effects. This approach is applied herein by providing a 
linear model at each system operating point. This will be dis- 
cussed next. 

Conveniently, an earlier parameter identification method 
for a linear system based on minimizing a cost function 
based on the error between a measured transfer function and 
a model transfer function which is described, for example, 
by Sanathanan and Koemer 33 is available. Miles 34 used this 
frequency domain type approach of minimizing a quadratic 
cost function based on the error between a measured boiler 
transfer function and a model boiler transfer function to 
model a system with a time delay and a similar type of 
approach is used on the system identification problem dis- 
cussed in this paper. 

In the next section the engine noise data are presented. 
Then the linear system theory will be discussed. Next, the 
acoustic spatial correlation model used herein to account for 
the enclosure will be presented. Then, results are presented 
and discussed. Finally, some conclusions are presented. 

II. ENGINE NOISE DATA AND ANALYSIS 

The engine test program and preliminary data analysis 
are described next. 

A. Engine test data 

The NASA/Honeywell static engine test program was 
conducted at Honeywell’s San Tan outdoor acoustic test fa- 
cility using a Honeywell TECH977 engine (Fig. 1) and the 
results are described in a report by Weir. 35 The dual-spool, 
turbofan engine has a direct drive, wide chord fan connected 
by a long shaft to the low-pressure turbine spool and a high- 
pressure compressor connected by a concentric short shaft to 
the turbine high-pressure spool. The fan diameter is 34.2 in. 
The combustor design is a straight-through-flow annular ge- 
ometry with 16 fuel nozzles and 2 igniters. Data obtained for 
one configuration in the test program are analyzed in this pa- 

Sensors in aft fan High-temperature sensors 

bypass access panels with air cooling at turbine exit 



with air cooling in 
combustor igniter port 

FIG. 1. (Color online) Honeywell TECH977 engine. 
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TABLE I. Spectral estimate parameters. 


Parameter 

Value 

Segment length (data points per segment), N 

32 768 

Sample rate, r s , samples/s 

65 536 

Segment length, Td = N/r s , s 

0.500 

Sampling interval. At — l/r s , s 

1/65 536 

Bandwidth resolution, B e — A f— 1 /Td= r s /N, Hz 

2.0 

Upper frequency limit, /„ = 1 /2At=r s /2, Hz 

32 768 

Number of independent samples, n s 

128 

Overlap 

0.50 

Sample length, r tota i, s 

«70 


per. The engine-internal instrumentation in this configuration 
included high-temperature pressure sensors with air cooling 
in a combustor igniter port identified herein as CIP1 (sensor 

9) and at the turbine exit sensors identified as T551 (sensor 

10) and T552 (sensor 11). Pressure time histories for the 
T551 turbine exit sensor are used herein. 

The data acquisition system had a sampling rate of 65 
536 Hz and a sampling duration of roughly 70 s. The spectra 
were calculated using a 50 percent overlap. This permitted 
data reduction using «254 overlapped ensemble averages at 
a bandwidth resolution of 2 Hz. Further signal estimation pa- 
rameters are shown in Table I. 

B. Combustor-turbine exit acoustics 

All the spectra and cross-spectra are estimated using 
Welch’s nonparametric method which is based on averaging 
multiple windowed periodograms using overlapping time 
sequences. 36 ’ 37 Using these spectra and cross-spectra, the 
magnitude squared coherence is calculated to measure the 
similarity of the amplitude variations at particular frequen- 
cies. The " accent will be used to denote the statistical basis 
of a variable. This is done to avoid confusion with calcula- 
tions of coherence using a single segment or block which 
yield a coherence of unity. 

The post-combustion (post-flame) residence time delay 
due to convection of entropy at the flow speed in the com- 
bustor identified by Miles 14 ' 15 also shows up in combustor 
sensor signal/turbine exit signal cross-spectrum phase angle 
plots. Based on the slope of the phase angle plot shown in 
Fig. 2, the time delay is slightly greater than 3.05 ms. In this 
paper, in order to study the enclosure sound field, the time 



FIG. 2. (Color online) Cross-spectrum phase angle between combustor and 
turbine exit signal in a turbofan engine at 48 percent maximum power (flight 
idle) using sensors CIP1(9) and T551(10) as measured (unwrapped). 
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delay is reduced by properly aligning the time histories 
by 200 time steps or about 3.05 ms (3.05176 ms 
= 200 x 1000/65 536). The alignment is used to reveal in 
more detail the interference pattern in the phase angle plots 
created by spatial correlation in the ambient sound field at 
low frequencies. The transfer function will be identified 
using the aligned signal cross-spectra. 

The combustor sensor signal/turbine exit signal magni- 
tude squared aligned coherence (MSC) function is shown in 
Fig. 3 for the 48 percent maximum power case. Plots using a 
linear and logarithmic MSC function scale are presented. 
The MSC function is given by 


Yx,y = 


I ,y I 

G , , G yj 


( 1 ) 


Also shown in Fig. 3(b) are coherence thresholds calcu- 
lated from 


i,w = l-a-p) 17 ^, (2) 

where this formula determines a E-percent threshold confi- 
dence interval using the number of independent data seg- 
ments/blocks, n s . used in the Welch’s periodogram method 
spectral estimator (see Miles 12 ’ 13 ). The 95 percent confi- 
dence interval based on n s =128 independent samples is 
0.0233. However, the spectra are calculated using a 50 per- 
cent overlap to reduce the variance and the 95 percent 
threshold confidence interval based on n s = 273 samples is 
0.0109. This horizontal line is also shown in Fig. 3(b). These 
indicators show the MSC function is reliable up to about 
400 Hz. However, as shown in Fig. 3(a) the MSC function is 
above 0.4 only in a region below 100 Hz. 



(a) Frequency, Hz 


— aligned 

— deliberately unalianed 



(b) Frequency, Hz 


FIG. 3. (Color online) Coherence using sensors CIP1(9) and T551(10) at 48 
percent maximum power (flight idle), (a) Linear scale, (b) logarithmic scale. 


Miles 12 ' 13 shows that instead of relying on the confidence 
interval given by Eq. (2), which is based on a statistical 
theory, to obtain a threshold value for y 2 n (« s ), one can use a 
deliberately unaligned time history to create the threshold 
value. If one of the time histories is shifted by a time delay 
more than the segment/block length, T d = N/r s , then the two 
time histories are totally independent unless tones are present. 
This deliberate decorrelation establishes a coherence thresh- 
old and also identifies any tones in the signals. Shifting the 
signals by this time delay removes the coherence of random 
noise but leaves the coherence of periodic functions which 
are sometimes identified as hidden periodicities, concealed 
spectral lines, or undamped sinusoids in noise. The deliber- 
ately unaligned coherence is also shown in Fig. 3(b). Note 
that the higher statistical confidence interval based on the 
number of independent records ( n d = 128) is a more conserva- 
tive estimate of the measured coherence threshold for the 
Honeywell TECH977 turbofan engine. Miles 13 previously 
found that with a 50 percent overlap the higher statistical con- 
fidence interval based on the number of independent records 
is a more conservative estimate of the measured coherence 
threshold for measurements made in the Pratt & Whitney 
PW4098 turbofan engine. Miles 1 1 showed that for a Gaussian 
distribution and using a 50 percent overlap to reduce the var- 
iance, a coherence threshold calculated with n a = 2 ns did 
work well. Consequently, simulations of internal turbofan 
engine noise with the assumption that it has a Gaussian distri- 
bution could lead to incorrect conclusions. The statistical co- 
herence threshold can be used with confidence since it can be 
compared with a measured coherence threshold. The coher- 
ence value is below the 95 percent statistical confidence inter- 
val above 400 Hz. Consequently, this is the upper frequency 
limit for which data are available for analysis using a linear 
system model. The linear system theory will be discussed 
next. Then, the enclosed ambient sound field theory will be 
developed using an approximate analytical formulation. 

III. UNENCLOSED LINEAR SYSTEM THEORY 

The methods used herein are based on system theory 
developed for linear systems with random inputs as 
discussed by Bendat and Piersol. -4 In this paper, the linear 
system theory discussed is in the frequency domain. The out- 
put spectral density function, G Vj3 ,, and the cross-spectra den- 
sity function, G xy , is related to an input spectral density 
function, G xx , through frequency response function, H xy (f), 
representing the turbine as 

Gy,y = \H x ,y(f)\ 2 G X)X (3) 

and 

G x ,y = l-h.y(f)G x , x , (4) 

where x is the input signal from the combustor pressure sen- 
sor CIP1 and y is the output signal from one of the turbine 
exit pressure sensors. 

The cross-spectral density and the frequency response func- 
tions are complex valued quantities, which can be expressed in 
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terms of a magnitude and an associated phase angle. This will 
be expressed herein using complex polar notation. 


GxM) = \Gx, y (f)\ exp 


-m) , 


(5) 


H x ,y(f) = \H Xt y(f)\exp[-j\l/(f)]. 


( 6 ) 


Before plotting the cross-spectral density phase angle, phase 
unwrapping is applied to the phase angle to avoid a jump of 
2n in the phase caused by the ATAN2 function. 17 

The system under consideration has a combustion noise 
input with a measured spectrum, Gg 9 (f), which includes 
acoustic and hydrodynamic components. The system meas- 
ured output quantities are assumed to be related as follows: 


<%, 10(f) = H %10 (f)^ 9 (f), 

G?o, 10 (f) = H 9,m (f) 2 G n 9,9(f) 


(7) 

( 8 ) 


and the measured MSC 


79,10 — 


J 9.10 


s^m f'm 

°9,9 ^10,10 


(9) 


where m indicates noise may be included in the measured 
quantities. Only measurements with the combustor sensor, 
CIP1 (9) and turbine exit sensor, T55 1 ( 1 0), are discussed 
herein. The unknown that will be identified is the turbine fre- 
quency response function, H™ w (f) at a range of operating 
conditions. 


IV. MODEL EQUATIONS 

The plant being modeled is the turbine and the sensor 
spatial correlation in the ambient noise field in the turbofan 
engine core. The input to the plant is the total pressure signal 
and the measurement made in the combustor is of the total 
pressure signal. Consequently, the available input auto- 
spectrum is that of the total pressure signal. A five parameter 
model is used for the transfer function. However, an eight 
parameter model is used for the enclosure sound field. Model 
parameters are adjusted so that agreement will be reached 
with the aligned measured cross-spectrum between the com- 
bustor and a turbine exit pressure sensor over the frequency 
range from 0-400 Hz and with the measured auto-spectra at 
the turbine exit sensor where the coherence is greater than 0.4 
from 25-100 Hz. In the region from 25-100 Hz the hydrody- 
namic noise is reduced and the turbine exit spectrum is not 
obscured by hydrodynamic noise. However, in this region it is 
obscured by the enclosure sound field. Normally, one might 
just fit the cross-spectra and compare the predicted and meas- 
ured auto-spectra at the turbine exit. However, the enclosure 
sound field shows up in the turbine exit auto spectra in the 
25-100 Hz range so that region is used in the curve fit. 

Next, the measured auto-spectra at the turbine exit and 
the estimated auto-spectra at the turbine exit will be com- 
pared and adjustments are made to the enclosure sound field 
model until a satisfactory result is obtained. The general idea 


is to show that the data indicate the transfer function inter- 
pretation is reasonable if the enclosure sound field is taken 
into account. To aid in physical interpretation, a standard 
template parametric model form will be used. The model 
form examined in this paper is in a parametric reduced order 
frequency domain representation. The parameters depend 
nonlinearly on the operating point. However, at each operat- 
ing condition the system will be assumed to be linear and the 
same parametric form will be used so that source separation 
will be obvious. The nonlinear operation is then described 
by a linear model at a range of observed operating points 
each identified by a set of parameters. This paper discusses 
measurements at the 48, 54, 60, 71, 87 percent maximum 
power settings and at a maximum power setting. The transfer 
function of the turbine is modeled using the simplest possi- 
ble dynamic system which is one which satisfies a first order, 
linear differential equation with a time delay for the indirect 
combustion noise, and no time delay for the direct combus- 
tion noise. Systems with time delays are in a class of systems 
discussed in books on linear time invariant system theory 
such as the text by Schetzen. 3x 

The objective in this section is to separate the combus- 
tion noise at the turbine exit from broadband noise from 
other sources and hydrodynamic pressure perturbations. The 
frequency region from about 0-50 Hz is a region where the 
measurements are masked by the ambient enclosure sound 
field which can be determined since the transducers are situ- 
ated within an acoustic wavelength of one another. An addi- 
tional ambient enclosure noise field model based on diffuse 
noise field theory is needed to match the measurements, and 
it will be discussed after the linear system theory sections. 

In this model any upstream propagating indirect com- 
bustion noise and any direct combustion noise reflected 
upstream from the turbine entrance are included in the com- 
bustion auto spectrum G xx (f). 

Consider an input signal x(t) with a spectrum G x Jf) for a 
system with transfer function H xy (f) and output signal y(t). 
Now that the statistical nature of the calculations is estab- 
lished, the" notation is dropped for simplicity. Then the cross- 
spectrum is given by Eq. (4). For the turbofan engine, the input 
to the turbine is the direct acoustic signal, xjt). and the time 
delayed entropy signal, x,(f), with a post-combustion residence 
time delay in the combustor of x pc . The time delayed signal 
may represent in addition to a temperature fluctuation moving 
with the flow any other disturbance moving with the flow such 
as a density or vorticity fluctuation. A turbine-combustor- 
tailpipe noise system diagram is shown in Fig. 4. Note the 
presence of a combustor entropy noise, N h combustor hydro- 
dynamic noise, N d , and turbine exit tailpipe noise, N T . 

The system model is constructed as follows. The output 
turbine noise signals are the direct and indirect combustion 
noise y d (t) and y,(t) as shown in Fig. 4. The output signal meas- 
ured at the turbine exit, y, is the sum of the direct combustion 
noise signal, y d , and the indirect combustion noise signal, y b 

7 = yi + yd- ( 10 ) 

The direct combustion noise cross-spectral density, G Xdyd , is 
a product of the direct combustion noise turbine transfer 
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Consequently, using Eqs. (13) and (14), the measured cross- 
spectral density at the turbine exit is given by 


Xi 


T 


yt 


Xd 


-G- 

N d l 


H i e iwT ' ,c 

1 

H d 

^Tailpipe '‘-t 

! 

y d 


FIG. 4. Turbine-combustor-tailpipe noise system. 


function, H d (f), and the direct combustion noise input 
spectral density, G Xd , and an acoustic travel time delay func- 
tion which is included for completeness 

G, dyd = H d (f)e' 2nflA G Xd . (11) 

When discussing auto spectra, the notation used is some- 
times simplified so that we write G Xd for G XdXd . 

The indirect combustion noise cross-spectral density, 
G Xdyn is a product of the indirect combustion noise turbine 
transfer function, //,(/), the time delay factor, e' lp , and the 
input indirect combustion noise spectral density, G x . , which 
corresponds to an equivalent fluctuating entropy spectral 
density function. 

G Xdy> = H i (f)e' <p G xn (12) 

where cp = (2nfx pc + (p o n/\80) I8O/71 is measured in 
degrees and t pc is the post-combustion residence time in the 
combustor. The pressure sensors are not phase matched and 
the phase angle (p a accounts for this. The two turbine exit 
sensors are assumed to be phase matched. 

The direct combustion noise and the entropy noise have 
the same origin in the combustion process. Consequently, it 
is assumed that at the turbine exit the resulting spectrum 
from the input entropy fluctuation spectral density can be 
modeled by an equivalent input acoustic spectral density and 
a transfer function 

Hi(f)e i<p G Xi = H(fy ,p G Xd . (13) 

This assumption is possible since the post-combustion time 
delay is observable in the cross-spectrum between the com- 
bustor and the turbine exit. Consequently, the observed indi- 
rect combustion noise and the direct combustion noise are 
coherent. 

The direct combustion noise turbine transfer function 
H d (f) in Eq. (11) is assumed to have a representation vH(f). 
Where v is a measure of the ratio of the direct combustion 
noise to the indirect combustion noise, 

H d (f) = uH{f). (14) 


G Xdiy = G Xdyd + G Xdy , = H d (f)^G Xd + H i (f)e' <l ’G Xi 
= H{f){^ + ve> 2 ^)G Xd . (15) 

In order to create a more desirable cross-spectmm phase angle 
variation with frequency so the enclosed sound field can be stud- 
ied, a time delay t„ is removed by properly aligning the time his- 
tories by 3.05 ms or 200 time steps (3.05176 = 200*1000/ 
65 536). The resulting cross-spectmm phase angle variation with 
frequency permits the determination of the sensor pair spatial 
correlation function with greater accuracy. 

The alignment time delay, t 0 , is known. This time delay 
is removed by a signal processing alignment procedure leav- 
ing the aligned cross-spectral density, U Xd<y , where 

IJ — e~ jlnfz °G 

u *i,y ‘ KJ *d,y 

= H(f) (e^- 27 ^) + ve i 2 nf(TA ~ z '‘ ) G Xd . (16) 


The combustor entropy noise, N h and combustor hydrody- 
namic noise, N d , are assumed to be independent of each 
other and independent of the tailpipe noise, N T . 


GN;N d — Gn^t — GnjNt ~ 0 . 


(17) 


The frequency response function of the turbine is mod- 
eled using the simplest possible dynamic system which 
is one which satisfies a first order, linear differential 
equation 


H(f) 


K 

1 +j2nfz T ’ 


(18) 


where K will be used to measure the indirect combustion 
noise created as a fraction of the noise in the combustor. 

This approach requires the input signal which is 
the measured combustor acoustic auto-spectrum, G99, 
be free of noise. However, the input signal is contami- 
nated by hydrodynamic noise and higher acoustic 
modes that are cutoff in the turbine. In addition, indi- 
rect combustion noise propagating upstream and direct 
combustion noise reflected from the turbine entrance 
may be present. The presence of noise limits the fre- 
quency range of the measured cross spectrum, G™ 10 to 
400 Hz because the cross-spectrum is not as affected 
by noise as long as the coherence is above the coher- 
ence threshold which is reached near 400 Hz as shown 
in Fig. 3. The frequency range of the measured turbine 
exit auto-spectrum, G”' 0 10 , unaffected by noise is lim- 
ited to 100 Hz based on the coherence function shown 
in Fig. 3 because the coherence falls below 0.4 at fre- 
quencies above 100 Hz. 

Since a frame of reference adjustment has been made 
using the alignment time delay, x,„ the estimated form of the 
free field system transfer function, H^ l0 (f), is 


'T/E, 

rl 9,10 — 


(eiiv-ZxfXo) _|_ 

1 + j2nf t t 


(19) 
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Hl w if)=KUl w , ( 20 ) 

where K is the coherent indirect combustion noise turbine 
transmission creation factor and Kv is the coherent direct 
combustion noise turbine transmission loss factor. 

The unenclosed transfer function factor 77 9 i0 will be 
used in the section discussing the enclosed transfer function. 
Consequently, it has a separate terminology. In the next sec- 
tion, the unenclosed transfer function Hg W (f) is discussed to 
specify relate the measured cross-spectra and auto-spectra to 
the proposed model if no enclosure was present. 


VI. SPATIAL CORRELATION COEFFICIENTS FOR 
AMBIENT CORE NOISE 

The enclosed space is shown in the cutaway engine dia- 
gram shown in Fig. 1. The dimensions are proprietary. Three 
elementary spatial correlation coefficients for ambient 
enclosed noise fields are discussed by Cook et al. 20 and 
Nelisse and Nicolas. 26 They are derived as follows by Cook 
et al. 20 except that in this derivation the assumption is made 
that at/=0 the sensors are not in phase. The spatial correla- 
tion between the sound pressures at two sensors A and B 
separated by a distance r m in a sound held is 


V. UNENCLOSED SOUND FIELD 


If free held data were available because the enclose was 
absent only a six parameter model would be used for the 


estimated transfer function to ht measured values of 


U[ 


9,10 


,|o' sncd , an( j Qm> ^ This would be done using estimated 
values as follows: 


101og(|£/f ilo |) = 101og 




1 +j2nf-z T 
+ 101og(/q + 101og(G" , 9 ), 


(21) 


T e \f) = tan” 1 10 (/)), K(//f I0 (/))} 180 /tt, 

(22) 


10 log(Gf 0 , 0 ) = 20 log 


( e /(<?>-2it/T D ) _|_ l , e -j2nf(z A -z 0 )) 


+ 20 log(A') 


1 + j2nfx T 

10 log(G 99 ), 


(23) 


where 


r/E 

u 9,10 


is the estimated free held combustor/turbine 


E aligned 
9,10 


is the esti- 


exit pressure cross-spectrum magnitude, r/j. 
mated free held combustor/turbine exit pressure cross- 
spectrum phase angle, and G h U) 10 is the estimated free held 
turbine exit pressure auto spectrum. The six parameters are K, 
v , x T , x pc , r A , and q> a . If the turbine exit noise is from the 
direct combustion noise source, vK is the turbine noise attenu- 
ation factor. However, if the turbine exit noise is from the 
indirect combustion noise source, K is an indirect combustion 
noise production/creation factor. The value of v indicates the 
contribution of direct combustion noise. The value of x T in 
addition to changing the phase angle also adjusts the slope of 
the magnitude of the estimated cross- spectral density function, 

r/E 

u 9,10 

The value of x A is much smaller than the value of t D 
since the speed of sound in the combustor greatly exceeds 
the how velocity. Consequently, t a will be neglected and 
only hve parameters deline the free held (missing enclosure) 
system. However, the measurements are masked by the am- 
bient enclosed sound held in the low frequency range which 
can be determined because the transducers are situated 
within an acoustic wavelength of one another in a closed 
space. The next section discusses the sensor pair spatial cor- 
relation model for ambient enclosed core noise. 


R = 


{pm) 

\A «>«» ’ 


(24) 


where pi and p 2 are the pressures at points A and B, 
respectively, and the angular brackets denote long time 
averages. 

If px and p 2 are both sinusoidal of the same frequency, 
but differ in phase by /i, then R = cos(/i) where 
pi = A cos (cat) and p 2 = B cos (cot — ft) since 

__j =s( ^ )=cos(w . (25) 


R = 3? 


Suppose a plane wave of wavelength X passes over the 
points A and B with i/r being the angle between the nor- 
mal to the wave fronts and the line AB. The phase dif- 
ference between the two points is = 2m "' c °^\ 
Consequently, the two points are in phase, i.e., fi = 2n if 
r, n = X and i/r = 0° and out of phase if t/r = 90°. Then the 
one dimensional ambient sound field spatial correlation 
for in phase sensors is = cos(/f) = cos(£r m cos(i/r)). 
However, the sensors used herein are not necessarily in 
phase at /= 0, so the one dimensional ambient core 
noise sound held spatial correlation used is Ri = cos (/f) 
= cos (kr m cos(i J/) — fl 0 ). 

If the sound held is random in two dimensions so that 
the normals to the wave fronts all lie in the same plane and 
all directions have equal weights, then the two dimensional 
ambient sound held spatial correlation is 


f2n 


R 2 = 


2n 


cos(Ar,„ cos(i/r) — , d 0 )dxlj = cos(-(? 0 )Jo {kr m ), 


(26) 


where ,/ 0 the Bessel function of zero order. 

If the sound held is random in three dimensions and all 
directions have equal weights, then the three dimensional 
isotropic ambient sound held spatial correlation is 


R 3 


1 

471 


f27I 


dcj) 


cos [kr m cos(i J/) 


cos(t? 0 ) 


sin(£/- m ) 

ki .I 


— d 0 ) sin(i/7)t/i// 


(27) 


The isotropic spatial correlation coefficient is used to evaluate 
the diffuse sound held in diverse spaces and reverberant spaces 
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such as rooms by Cook et al. 20 and Nellisse and Nicolas, ° 6 and 
wind tunnels by Piersol ’ and Bendat and Piersol." 

In an isotropic noise field the directional density of the 
noise is uniform so that the same amount of noise arrives 
from all directions. The core noise ambient field spatial cor- 
relation model will be constructed using a nonisotropic two 
dimensional noise field spatial correlation function. In addi- 
tion, an anisotropic one dimensional noise field will be used. 
This is similar to the use of an isotropic three dimensional 
noise field and an anisotropic surface noise source field in 
the study of ambient sea noise by Cron and Sherman." 1 

The models will be specified by an a ambient core noise 
travel time t using the relationships 


<m gned (/) = tan- 1 {s(^ 10 (/)),3?« 10 (/))}l80/7t, 

(33) 

10 log(Gf 010 ) = 20 log(|w£ 10 |) + 20 logOST) 

+ 10 log(G^), (34) 


where the enclosure sound field is a fraction of the combus- 
tion sound field, G99, and K is an appropriate attenuation 
factor. 


VIII. PARAMETER EVALUATION PROCEDURE 


c 

c = If, 


kr m = '■ 


2nr r 


2nfr K 


= 2 nfx, 


since while distances and velocities are unknown, the fre- 
quency is known and the appropriate time constant t can be 
estimated from the data. This approach avoids limitations 
imposed by the lack of proprietary geometry and flow 
information. 

In this paper, we define a one dimensional anisotropic 
core spatial correlation function where we assume the nor- 
mal of the wave fronts is parallel to the line between the sen- 
sors and consequently t/r = 0. Three parameters are used 
since the sensors may not be in phase. Two of these spatial 
correlation functions are used herein 


Ci (/') = «ci cos(27t/T cl - t? c i), (28) 

C2 (/') = «c2 cos(27t/i c2 - fid), (29) 

and a two dimensional ambient turbine exit spatial correla- 
tion function where we use two parameters 

'/(/) = a T EJo{2Tift TE ). (30) 

The total spatial correlation, p(f) is 

p(f) = *?(/) + Ci (f) + C2 if)- 


VII. ENCLOSED SOUND FIELD 

This model is designed to fit the large scale structure in 
the interference pattern in the phase angle due to the entropy 
waves which have longer time delays. 

We then have for the ambient acoustic closed space core 
noise transfer function, cross-spectrum, and turbine exit 
auto-spectrum, which are compared to the measurements 


^ 1 0 = ^,i 0 + C 1 (/) + C 2 (/)+t/(Z), (31) 


10 



= 10 1og(w£ 10 ) 

+ 10 log(tf) + 10 log(G^), 


(32) 


Model parameters were first selected and adjusted so 
that the auto-spectra and cross-spectra created by the system 
model were in fair agreement with the measured values. 
Next, an appropriate cost function is used to fine tune these 
model parameters. The cost function is outlined in the 
Appendix. The cost function uses factors to place more 
weight on the least square error values which are considered 
more reliable. A higher reliability is placed on the measured 
phase angle from 0-250 Hz, the measured cross-spectrum 
magnitude, and the turbine exit auto-spectrum from 0-100 
Hz than the measured phase angle from 252—400 Hz. A solu- 
tion method that provides an optimum solution in a least 
squares sense without derivatives is used. The search tech- 
nique used in this study is described by Powell 39 and a FOR- 
TRAN computer code for this algorithm is given by Shapiro 
and Goldstein, 40 and Kuester and Mize. 41 The model param- 
eters determined for the 48, 54, 60, 71, 87 percent of maxi- 
mum power operating condition cases and one maximum 
power case are shown in Table II, and selected parameters 
are shown in Fig. 5. The exact power setting of the maxi- 
mum power case is unknown. To make it possible to plot the 
maximum power parameters, the maximum power case is 
arbitrarily assigned a value of 99 percent of maximum 
power. 

Figure 5 shows parameters related to the turbine and 
flow system. Figure 5(a) shows the ratio of direct to indirect 
coherent combustion noise v is about 26 percent at the low- 
est power setting and then drops from 12 percent at 50 per- 
cent of maximum power to 5 percent at the maximum 
power setting. Figure 5(b) shows the post-combustion resi- 
dence time in the combustor x pc ranges from about 3.4 ms 
at low engine speed to about 2.9 ms at the maximum power 
setting. 

Examining Table II shows the turbine acoustic attenua- 
tion/creation level, 10 log (K), ranges from —5 to —9 dB. 
This means the loss/production of acoustic power across the 
turbine that can be attributed to the transfer function ranges 
from 10 dB to 18 dB. The level of the coherent noise at the 
turbine exit attributed to the transfer function is 10 dB to 18 
dB below the level of the combustion noise. An additional 
loss is apparent which is attributed to phase cancelation due 
to the spatial correlation in the ambient core noise field. This 
can be more than 6 dB. 

Examining Table II shows the system behaves like a 
first order linear system with a time constant, x T , ranging 
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TABLE II. Model parameters. 


Parameter 
Percent maximum 

Value 

48 

54 

60 

71 

87 

99 

jAg) 

0.2585567 

0.09519884 

0.1227427 

0.08256044 

0.08220044 

0.05111235 

to logWfe) 

-8.38677 

-5.007134 

-7.788622 

-9.343841 

-7.906158 

-7.422281 

x T (€ 3 ), s 

0.001595392 

0.0006247073 

0.001532724 

0.001444117 

—0.0009414114 

— 3.963765e-09 

V (£ 4 ), ms 

3.437903 

3.446463 

3.293785 

3.047519 

3.172503 

2.940178 

q> 0 (e 5 ), deg 

0.2492384 

-2.887404 

-1.424596 

-6.793756 

-4.927228 

-18.28111 

a-TE (^6) 

-0.9847035 

-0.3625233 

-0.5555566 

-0.2957081 

-0.5565302 

-0.6333497 

x TE (£ 7 ), ms 

7.659317 

7.83757 

7.748561 

7.443147 

7.839088 

8.172998 

1 (eg) 

-0.2224676 

-0.5898645 

-0.3994165 

-0.3836406 

-0.4160725 

-0.3804713 

x c 1 (e 9 ), ms 

6.013163 

2.565741 

3.541956 

2.361771 

3.108431 

3.203457 

tfc, (eio), deg 

89.12107 

-0.2313568 

29.9959 

360.0005 

61.16767 

90.33118 

«c2 (Gl) 

-0.3837000 

-0.9980998 

-0.3701099 

-1.999444 

-0.4794316 

-0.8302187 

x c i (ei 2 ), ms 

3.986773 

1.334677 

2.658769 

0.1102978 

1.634171 

1.037572 

$C 2 , deg 

90.68719 

88.36214 

90.15674 

274.7305 

90.09822 

86.5299 


from 0 to 1.5 ms. The time constants for the one dimensional 
ambient core noise spatial correlation functions, 4 , if) and 
£ 2 (/), which are t,. and z C2 , and the two dimensional ambient 
core noise spatial correlation function, 11 (f), which is x TE are 
shown in Table II. At low power settings the three dimen- 
sional ambient core noise spatial correlation function time 
constant x TE is about 8 ms. At low power settings the one 
dimensional ambient core noise spatial correlation function 
time constant r Cl and r C2 ranges from 0 . 1-6 ms. 

IX. RESULTS 

Results for the 48, 54, 71, and 87 percent maximum 
power case are shown in Figs. 6-9. Two figures are used to 
compare the estimated acoustic cross-spectrum magnitude 
and phase with measured values over the frequency range 
from CM-00 Hz. These figures also show the combustion 




FIG. 5. Turbine/flow system parameters as a function of percent maximum 
power, (a) Ratio of direct to indirect combustion noise, (b) post-combustion 
residence time in combustor. 
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region auto-spectra and the measured and estimated turbine 
exit auto-spectra. 

The large scale variation of the cross-spectrum phase 
angle with frequency is captured by the ambient core noise 
spatial coherence function used herein as shown in (a) of 
Figs. 6-9. 

The estimated and measured cross-spectrum magnitude 
and phase angle shown in Figs. 6-9 are in good agreement 
from 0-400 Hz. The estimated and measured turbine exit 
auto-spectrum shown in (b) of Figs. 6-9 are in fair agree- 
ment from 0-100 Hz in these figures. Above 100 Hz, the 
auto-spectra results at the turbine exit are obscured by hydro- 
dynamic and uncorrelated acoustic noise. 

The spatial correlation function, p(f), and its compo- 
nents [ 11 (f), (i(f), ( 2 (/)] which are used to model the data for 
each case are shown in Figs. 10 and 11. For clarity, since at 



FIG. 6 . (Color online) Cross-spectrum turbofan engine at 48 percent maxi- 
mum power (flight idle) using sensors CIP1(9) and T551(10) from 0^100 
Hz. Phase angle shown using alignment of 200 counts or 3.05 ms. (a) Phase 
angle, (b) magnitude. 
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FIG. 7. (Color online) Cross-spectrum turbofan engine at 54 percent maxi- 
mum power (approach) using sensors CIP1(9) and T551(10) from 0-200 
Hz. Phase angle shown using alignment of 200 counts or 3.05 ms. (a) Phase 
angle, (b) magnitude. 

/= 0 , the components and the spatial correlation function are 
negative, each item has been multiplied by —1 before plot- 
ting. The observed 9 dB to 10 dB attenuation of the combus- 
tion noise shown in plots of the cross-spectrum magnitude in 
Figs. 6 (b) — 9(b) is normally attributed to the turbine. Exami- 
nation of Figs. 10 and 11 shows the actual turbine attenua- 
tion can be less by 2 to 6 dB or more at some frequencies. A 
good part of the attenuation measured may be due to destruc- 
tive interference due to the spatial correlation in the enclosed 
core noise sound field of the turbofan engine. As a result, the 



FIG. 8. (Color online) Cross-spectrum turbofan engine at 71 percent maxi- 
mum power (cutback) using sensors CIP1(9) and T551(10) from 0—400 Hz. 
Phase angle shown using alignment of 200 counts or 3.05 ms. (a) Phase 
angle, (b) magnitude. 


FIG. 9. (Color online) Cross-spectrum turbofan engine at 87 percent maxi- 
mum power (takeoff) using sensors CIP1(9) and T551(10) from 0^4-00 Hz. 
Phase angle shown using alignment of 200 counts or 3.05 ms. (a) Phase 
angle, (b) magnitude. 

transfer function attenuation 10 log (K) shown in Table II 
ranges from —5 to — 9 dB. 

The total cost function value, C, and the components S), 
S i a , S 2 , and S 3 are defined in the Appendix and tabulated for 
each power setting in Table III. The components are least 
square fit values. The value of the least square tit error to the 
cross-spectrum phase angle from 0-250 Hz, Si, increases as 
the engine speed increases indicating the simple model used 
matches the measurements better at low engine speeds. This 
is also reflected in the total cost function C also shown in 



FIG. 10. (Color online) Spatial correlation function at 48 and 54 percent 
maximum power, (a) 48 percent maximum power (flight idle), (b) 54 percent 
maximum power (approach). 
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FIG. 11. (Color online) Spatial correlation function at 71 and 87 percent 
maximum power, (a) 71 percent maximum power (cutback), (b) 87 percent 
maximum power (takeoff). 

Table III. The value of the least square fit error to the cross- 
spectrum phase angle from 252^400 Hz, Si a , decreases as 
the engine speed increases indicating the simple model used 
matches these measurements better at high engine speeds. 
The least square fit error to the cross-spectrum magnitude, 
S 2 , does not vary much with engine speed. The least square 
fit error to the turbine exit auto spectrum, S 3 , is larger at the 
higher engine speed. 

X. DISCUSSION 
A. Estimated transfer function 

The most significant result is the good agreement 
between the aligned measured and estimated cross-spectra 
shown in Figs. 6-9. That this is achieved with only a five pa- 
rameter curve fit using as a transfer function of the turbine 
the simplest possible dynamic system which is one which 
satisfies a first order, linear differential equation with a time 
delay for the indirect combustion noise and with no time 
delay for the direct combustion noise is remarkable. The lin- 
ear system model produces good agreement between the 
measured and estimated cross-spectral density function from 
CM-00 Hz. In addition, the linear system model produces fair 


agreement between the measured and estimated turbine exit 
auto-spectral density function from 0-100 Hz where the sys- 
tem acoustic random noise is negligible. 

At some frequency in the 100—400 Hz frequency range 
the measured turbine exit auto spectrum and the measured 
combustor/turbine exit cross-spectrum magnitude approach 
the same level and there is no way to make judgments about 
the transfer function approach to predict the turbine exit 
auto-spectrum. Hydrodynamic or uncorrelated noise makes 
it impossible to compare the estimated turbine exit auto- 
spectrum in the 100 to 400 Hz frequency range. However, 
the results suggest that indirect combustion noise is most im- 
portant in the 0-100 Hz frequency range. The results show 
the importance of treating the turbofan noise source as a sys- 
tem rather than a collection of separate noise generating 
parts. The results provide insight into the part indirect com- 
bustion noise plays as a source low frequency core noise for 
this engine. 

The observed attenuation of the combustion noise 
shown is normally attributed to the turbine. The actual tur- 
bine attenuation is less by 2 to 6 dB or more at some fre- 
quencies. A large part of the attenuation measured is from 
destructive interference due to the spatial correlation in the 
core noise sound field of the turbofan engine. This underly- 
ing physics is not included in current combustion noise pre- 
diction schemes. The current combustion noise prediction 
capability applies over a very narrow design space. Conse- 
quently, new turbofan engine core configurations might be 
outside the current combustion noise prediction capability. 

B. Spatial correlation, interference patterns, and 
nonlinearity 

The correlation structure of the enclosed ambient core 
noise shown in Figs. 10 and 11 differs from that of a diffuse 
sound field which might be found in measurements in rever- 
berant rooms or wind tunnels. A diffuse sound field is one in 
which the time average of the mean-square sound pressure is 
the same everywhere and the flow of acoustic energy in all 
directions is equally probable. The model is anisotropic 
because the ambient core noise sound field is based on a two 
dimensional spatial correlation function and a one dimen- 
sional spatial correlation function with a sound field parallel 
to the microphone axis. An anisotropic formulation of spatial 
coherence is also used in the study of ocean noise spatial co- 
herence. For example, the ambient noise field in the ocean is 
analyzed as a superposition of an anisotropic noise field due 
to the ocean surface in an isotropic noise field by Cron and 


TABLE III. Cost function analysis, C = (Si * 1000 + S,„ * 200 + S 2 * 800 + S 3 * 1000)/201 . 


Cost 

Percent maximum 

Value 

48 

54 

60 

71 

87 

99 

c 

0.6250858 

0.5526794 

0.6700567 

0.8762981 

1.01022 

0.9716793 

Si 

0.03276443 

0.03940959 

0.04935671 

0.0527839 

0.05993345 

0.05132201 

Sla 

0.05608073 

0.04473534 

0.0387119 

0.0335156 

0.01195643 

0.01648785 

s 2 

0.07656069 

0.05218781 

0.07008126 

0.09860328 

0.1180152 

0.1209408 

s 3 

0.02037054 

0.02098165 

0.02151524 

0.03776569 

0.04631416 

0.04392558 
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Sherman. 21 These measurements also differ from those 
made in reverberant rooms or wind tunnels since the sound 
transducers are not phase matched. 

An alternative derivation of the one, two, and three 
dimensional spatial correlation relationships was made by 
Morrow 22 using the actual mode shapes of a rectangular 
room and some restrictive assumptions. He suggests that 
large deviations from the formula are possible if his assump- 
tions are not met. Because one assumption is that the modal 
density is high and this is not true at low frequencies, the 
deviations observed using these standard spatial correlation 
models are not unexpected. 

The results presented herein are not primarily an exam- 
ple of nonlinearities appearing in turbofan engine low fre- 
quency core noise. Instead, they are explained as being 
mainly due to an interference pattern caused by the spatial 
coherence of the sensor pair in the ambient core noise field. 
The interference pattern caused by the spatial coherence of a 
sensor pair in a diffuse noise field has been treated as part of 
reverberant room studies.- ’ - The literature on the spatial 

coherence of a sensor pair in the presence of turbofan core 
noise has a lacuna. 

However, the results in the 25-50 Hz region where the 
peak sound occurs are in better agreement at the low power 
settings and in poor agreement at the high power settings. 
This might be attributed to some combination of nonlinear 
effects and spatial coherence effects not identified with the 
simple model discussed herein. However, the post- 
combustion residence time measured herein seems to indi- 
cate the source of the sound in the low frequencies is due to 
indirect combustion noise. 

XI. CONCLUDING REMARKS 

The presence of coherent indirect and direct combustion 
noise in the core of this turbofan engine is observed. A low 
order model for the frequency response of the turbine to 
convected entropy and other disturbances moving at the flow 
velocity and direct combustion noise is successful in calcu- 
lating the combustor/turbine exit pressure cross-spectral 
density from 0-400 Hz and the turbine exit pressure auto- 
spectrum from 0-100 Hz at all speeds. The ratio of direct 
combustion noise to indirect combustion noise is about 26 
percent at the lowest power setting and then drops from 12 
percent at 50 percent of maximum power to 5 percent at the 
maximum power setting. 

An important result of the analysis of the turbine trans- 
fer function was development of a scheme to measure the 
post-combustion residence time in the combustor of a turbo- 
fan engine which might be useful in predicting or reducing 
thermal NO, and soot. Also, a measurement of the post- 
combustion residence time in the combustor might be useful 
in measuring the relative health of a gas turbine engine or 
doing real time performance monitoring. 

In order, to relate the transfer function model to the 
measurements studied, a spatial correlation function is 
developed for the enclosed ambient core noise field of the 
turbofan engine. Because proprietary geometry and flow in- 
formation is not available, the approach used is based on an 


approximate analytical formulation. In addition, in this pa- 
per in order to study the enclosure sound field, the convec- 
tive time delay observed in the combustor/turbine exit 
cross-spectrum phase is reduced by properly aligning the 
time histories. The cross-spectrum phase angle calculated 
after alignment is used because it reveals in more detail the 
interference pattern in the phase angle plots created by spa- 
tial correlation in the ambient sound field at low frequen- 
cies. The large scale structure interference pattern observed 
in the cross-spectrum phase angle measurements was repli- 
cated using this model. In addition, the model helps explain 
the measured noise reduction at low frequencies. The 
model shows this reduction is due to phase cancellation. 
The interference pattern can be observed in the cross- 
spectrum phase measurement. An anisotropic pair of spatial 
correlation functions is used to model the ambient core 
noise spatial correlation function. The observed attenuation 
of the combustion noise shown is normally attributed to the 
turbine. The actual turbine attenuation is less by 2 to 6 dB 
or more at some frequencies because a large part of the 
attenuation measured is from destructive interference due 
to spatial correlation in the core noise sound field of the tur- 
bofan engine. The advance demonstrated here is that by 
providing an enclosure sound field based on spatial correla- 
tion in the engine core using an approximate analytical for- 
mulation the resulting model is in improved agreement 
with measured acoustic cross-spectra phase. As a result the 
acoustic linear system transfer function model between the 
combustion noise and the turbine exit noise was to some 
extent validated. 
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APPENDIX: COST FUNCTION 

The cost function used is 

C = (Si* 1000 + S la * 200 + S 2 * 800 
+ Si * 1000) /201 , 


where 


/— 250 Hz 

E 

/= 0 


(f) 


1000 


/=400 Hz 
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/=100 Hz 


^3= E 

f=0 


log( 

G(f)u 3,10 

)-log( 

G(f)? 0 , 10 ) 

1 2 

10 



and the phase angles are unwrapped. 

The solution space is restricted by additional cost func- 
tions as follows: 


Cd(/' = o) < o, 

Cdif = o) < o, 

£i : 0 < v < 1.0, 
e 2 : 10 log(K) < -5, 
£3 : 0 < Tr < 0.01, 

£4 • 0 < tpc < 5, 

£5 : 0 < \cp a \ < 360, 
£ 6 : -2 < |at £ | < 0, 

£7 : xje unrestricted, 
e 8 : —2 < |a c i | < 0, 
£9 : t c i unrestricted, 
£to : Ki | < 360.0, 

£11 : —2 < |a c2 | < 0, 

£12 : t C 2 unrestricted, 
£ 13 : |t? c2 |< 360.0. 
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